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1. Introduction 


As a possible approach to the functional integral with complex measure, such as the one 


the complex Langevin method 


has 


encountered in the finite density QCD []]. 
attracted much attention in recent years. This recent interest was triggered mainly by the 
discovery of sufficient conditions for the convergence of the method to a correct answer [g, [ 7 ]. 
Reference [g] is a review on recent developments. Roughly speaking, if the probability dis¬ 
tribution of configurations generated by the Langevin dynamics damps sufficiently fast at 
infinity of configuration space, the statistical average over the configurations is shown to be 
identical to the integration over the original complex measure. R has been observed that, in 
systems for which the complex Langevin (CL) method converges to a wrong answer (such as 
the three-dimensional XY model |9|), this requirement of a sufficiently localized distribution 
is broken, typically in “imaginary directions” in configuration space. 

After the above understanding, a prescription in lattice gauge theory that makes the prob¬ 
ability distribution well localized was proposed in Ref. [To^; the prescription is termed “gauge 
cooling” and it proceeds as follows: The link variables in lattice gauge theory are originally 
elements of the compact gauge group SU(N). When the (effective) action is complex, how¬ 
ever, the corresponding Langevin evolution drives link variables into imaginary directions 
and link variables become elements of SL(N, C), a noncompact gauge group@ This evolu¬ 
tion tends to make the distribution wide in noncompact directions; in terms of the SU(N) 
Lie algebra, those noncompact directions are parametrized by imaginary coordinates. At 
this point, one notes that the definition of a physical observable that is invariant under 
the original SU(N) gauge transformations can always be tailored so that it is invariant 
also under the noncompact SL(N,C ) gauge transformations. The idea of the gauge cooling 
is that by applying the SL(N, C) gauge transformations appropriately along the complex 
Langevin evolution, one squeezes the distribution well localized so that the prerequisite of 


the convergence theorem [6j, [7]] is fulfilled without changing physical observables. 

In a one-dimensional gauge model and in the four-dimensional QCD with heavy quarks, 
it has been confirmed that the gauge cooling makes the distribution well localized and the 
complex Langevin method gives rise to correct answers [ltl]. More recently, this method 
was applied to the full QCD at finite density 11]. See also Refs. 12- 14J. One should note, 
however, that the Langevin dynamics itself is defined on gauge noninvariant variables (i.e., 
link variables) and also that the gauge cooling step cannot be regarded as a Langevin evo¬ 
lution that is induced by a holomorphic action; the latter is assumed in the convergence 


theorem [g, |7|. Strictly speaking, therefore, the convergence theorem does not apply when 
the gauge cooling is employed. The method should still be carefully examined in various 
possible ways. 

In the present paper, we apply the complex Langevin method in conjunction with the 
gauge coo ling to the two-dimensional lattice Yang-Mills theory which can be analytically 
solved [15-1171]. By doing this, we examine the validity of the method. The partition function 


i 


We will shortly describe the Langevin evolution of link variables. 
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of the two-dimensional Yang-Mills theory on the lattice is given by 


Z = 


/ 


Jd U Xtll 



( 1 . 1 ) 


where U X , L are link variables defined on a two-dimensional rectangular lattice, S is the lattice 
actionH 

S = E Tl ' [^01 (s) + U m (x)- 1 ] , (1.2) 

X 

and the plaquette variable is defined by 

Ufxvix) = U x ^U x ^jx t uU x+illJ U xl ,. ( 1 - 3 ) 

For simplicity, we assume that the gauge group is 517(2), that is, U X) ^ € SU{2) in the 
original integral CL3D- On the other hand, when the gauge coupling /3 is complex, the corre¬ 
sponding Langevin equation [Eq. (12.11) below] evolves link variables as elements of SL( 2, C). 
Thus, the distinction between U x ,n and becomes very important in the complex 

Langevin dynamics. For the convergence theorem in Refs, p, [?] to apply, the action 5 
that generates the drift force in the Langevin equation and physical observables must be a 
holomorphic function of dynamical variables; our above definitions Cd Q that entirely 
use U~^ not U x iAJ are chosen by this criterion. Note also that the plaquette action (11.21) is 
invariant under the 5L(2, C) lattice gauge transformations [such as the one in Eq. (12.41) ]. 

We consider the expectation value of the plaquette variable: 


<Tr[[/ 0 i (*)]> = 


J dU X)/l 


Lx,ii 


-,-S 


TV[C/ 0 i(x)]. 


(1.4) 


Even if the gauge coupling /3 is complex, this can be exactly computed by the character 


expansion 


15Ml7l]. Under periodic boundary conditions, one yields 


N d 


(Tt[(/ M (x)]) = - F —ln2, 2 = ^ 


n=1 


1 V 




(1.5) 


where I n (x) denotes the modified Bessel function of the first kind and V is the number of 
lattice points. 


2. Complex Langevin method and the gauge cooling 

The following procedures are basically identical to the ones adopted for the four-dimensional 
lattice QCD in Ref. 11] for example, although our two-dimensional pure-gauge system is 
much simpler. 


2 Throughout this paper, N = 2. 
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For the link variable, the Langevin equation with a discretized Langevin time t with the 
time step e is defined by 


U X M + e) = exp 


^ ^ ^ Aq (t) x 


UxAt)> 


( 2 . 1 ) 


where A a (a = 1, 2, 3) are Pauli matrices, r] atX ^{t) are Gaussian real random numbers of the 
variant 

)) — j (2-2) 

and D a x ^S is the drift force generated by the action S in Eq. OD; the derivative with 
respect to the link variable is given by 


Da,x,yf(U) = dsfie^U^) 


€=0 


(2.3) 


When the gauge coupling /3 is complex, the drift force becomes complex and the Langevin 
evolution evolves link variables as elements of SL( 2,C). 

The above complex Langevin dynamics tends to make the probability distribution function 
of link variables wide in noncompact directions of SL( 2, €). To squeeze the distribution well 
localized without changing gauge invariant quantities, we apply the following SL( 2, C) gauge 
transformation (this step is the gauge cooling) 


U, 


X 


K, = VxU x ,,v - + 1 A , 


where 


V x =e-"-'- K , /a — 2 TV 




and a > 0 is a real parameter. The distance defined by [18] 

X,(l 


(2.4) 

(2.5) 

( 2 . 6 ) 


measures how a SL( 2, C) gauge field is far away from the subspace of SU( 2) gauge fields. It 
is then straightforward to see that for a sufficiently small e, the gauge cooling (12.41) decreases 
or does not change the distance d. Note that in Eq. (12.51) is not a holomorphic function of 
link variables and thus the step (12.41) cannot be regarded as a part of the complex Langevin 
dynamics in which the drift force is generated by a holomorphic action; this fact prevents 


us from applying the convergence theorem 


to the above procedures. 


3. Result of numerical simulations 

We numerically solved the Langevin equation with the discretized Langevin time, Eq. (El]), 
on a V = 4 2 lattice. Periodic boundary conditions are imposed. The maximal size of the time 
step e we adopted was 0.001 and, when the drift force becomes large, we further reduced e 
“adaptively” according to the prescription in Ref. [191 ]. 

As the parameter a in the gauge cooling (12.51) , we tried both a = 1 and the adaptive choice 
(see Ref. 


Qad = A. o = ^Ei/«i + l 


D 


(3.1) 


Our numerical results did not show any notable difference in these two choices and we present 
the results with the latter choice in what follows. 
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To determine an appropriate rate of the gauge cooling (12.41) along the Langevin evolution, 
we observed the time evolution of the distance d (12.61) by changing the number of the gauge 
cooling steps per one Langevin update (12.11) . In Fig. [TJ for /3 = 0.4 + 2.0 i, we plotted the 
evolution of the distance d as the function of the Langevin time t by changing the number of 
the gauge cooling steps per one Langevin update as, 10, 30, and 100@ Since this plot shows 
the evolution including the Langevin stochastic dynamics, the distance d does not necessarily 
decrease. Since we do not see much difference for those three choices, we adopted ten gauge 
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Fig. 1 Evolution of the distance d (|2.6D for /3 = 0.4 + 2.0 i with various numbers of the 
gauge cooling steps per one Langevin update (12.lj) , 10, 30, and 100. 


cooling steps per one Langevin update. With this choice, the evolution of the distance d 
for various complex gauge couplings, f3 = 0.4 + 0.4 i, 0.4 + 2.0 i, 2.0 + 0.4 i, and 2.0 + 2.0 i, 
looks as depicted in Fig. [2j It appears that the gauge cooling is working perfectly for those 
complex gauge couplings, suppressing the evolution to noncompact imaginary directions. 

Now, we turn to the computation of the expectation value of the plaquette, Eq. (HD, by 
the complex Langevin method. Starting from a configuration of random SU(2) matrices, we 
discarded configurations until the Langevin time t = 11 for thermalization. Then 1000 con¬ 
figurations separated by At = 1 from t = 11 to t = 1010 are used to compute the expectation 
value. For typical values of the complex gauge coupling, we confirmed that the plaquette 
values between configurations separated by At = 1 practically have no autocorrelation. In 
Figs.[3HH we plotted the real and imaginary parts of the expectation value (ll.4|) obtained by 
the CL method. The error bars are statistical ones. The horizontal axis is the the complex 
phase 9 of the gauge coupling with the modulus 1.5: 

/3 = 1.5e i9 , 0 < 0 < 7t/2. (3.2) 


3 We define the Langevin time t such that it does not elapse during the gauge cooling steps. 

4 When the gauge group is SU( 2), the partition function (11.11) is invariant under /3 —ft. Because 
of this invariance and the complex conjugation, it is sufficient to consider the range of the phase, 

0 < 9 < tt/2. 
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Fig. 2 Evolution of the distance d (12.6p for various complex gauge couplings, f3 = 0.4 + 
0.4 i, 0.4 + 2.0 i, 2.0 + 0.4 z, and 2.0 + 2.0 i. The number of the gauge cooling steps per one 
Langevin update (12.11) is 10. 


The solid line curves are exact values given by Eq. (USD. We see that the complex Langevin 



Fig. 3 Real part of the expectation value (da obtained by the CL method and the exact 
value given by Eq. (11.51) . The horizontal axis is the complex phase 6 of the gauge coupling 
in Eq. (13.21) . 


method reproduces the real part fairly well, while it clearly fails to converge to the correct 
value of the imaginary part when the complex phase of the gauge coupling is large. 
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Fig. 4 Imaginary part of the expectation value (|1.4j) obtained by the CL method and 
the exact value given by Eq. (11.51) . The horizontal axis is the complex phase 9 of the gauge 
coupling in Eq. (HH- 


The gradation plot in Fig. [7] shows the relative error 

|(Tr[C„ I (x)]) CL -(Tr[C„ 1 (x)]) e ^| 


|<IV[V„i( I )]> m J (3 ' 3) 

on the first quadrant of the complex /3 planeH The (quadrant) circle in the figure is Eq. (13.21) 
along which Figs.EHUare plotted. Clearly, the relative error of the complex Langevin method 
becomes large when the complex phase of the gauge coupling becomes large. Four black 
crosses in the figure indicate complex gauge couplings we used in Fig. [2} the behavior in Fig. [2] 
thus suggests that the gauge cooling is correctly working for the region of the complex gauge 
coupling shown in Fig. 0 Nevertheless, the complex Langevin method shows large deviation 
from the correct value as in Fig. [H This is the main result of the present paper. 

A similar failure of the complex Langevin method for large complex j3 has been observed; 
see Fig. 8 of Ref. [201 ]. This result of Ref. [20( is, however, for a one-dimensional integral (not 
a gauge theory) and the validity of the gauge cooling, which is our main issue in this paper, 


is not relevant to this result of Ref. [20|. 


It is of interest how configurations generated by the Langevin dynamics distribute in 
configuration space. To give some idea on this point, in Figs. [6] [71 we present scatter plots 
of the plaquette variable (averaged over the lattice volume) for each configuration. Both 
cases, with and without the gauge cooling, are shown. Figure [6] is for j3 = 1.5 e* (°- 37r / 2 ) (i.e., 
9 = 0.37t/2) and corresponds to points in Figs. SHU with a relatively small complex phase. 
Figure [7] is, on the other hand, for /3 = 1.5 e* (°- 77r / 2 ) (i.e., 9 = 0.77t/2) and corresponds to 
points in Figs. [3}[I]with a large complex phase and with large deviation. Although there is 
a tendency when the complex phase of the gauge coupling is large for the distribution to 
become somewhat wider even after the gauge cooling, it is not clear from the scatter plot 


5 The block around the origin /3 = 0 is omitted from the plot because (Tr[?7oi(a:)]) ~ 0 for ft ~ 0. 
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Fig. 5 Gradation plot of the relative error (13.31) on the complex /3 (the gauge coupling) 
plane. The block around the origin /3 = 0 is omitted from the plot. The quadrant is Eq. (|3.2I) 
along which Figs. [3] HI are plotted. Four black crosses indicate complex gauge couplings we 
used in Fig. [2) 
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Fig. 6 Distribution of the plaquette variable averaged over the lattice volume. /3 = 
1.5e*(°- 37r / 2 ) (Q = 0.3tt/2). 


in Fig. [1] alone whether the distribution is so poorly localized as to break the prerequisite of 
the convergence theorem [g, 7*]. More detailed study is needed on this point. 

4. Conclusion 

In the present paper, we applied the complex Langevin method in conjunction with the 
gauge cooling to the two-dimensional lattice SU(2 ) Yang-Mills theory. Our intention was 
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Fig. 7 Distribution of the plaquette variable averaged over the lattice volume. /3 = 
1.5e*(°-W 2 ) (Q = 0.7tt/2). 


to examine the validity of the method by using this analytically solvable model. Somewhat 
unexpectedly, as shown in Figs.HHSJ we obtained strong numerical evidence that the method 
fails to converge to the correct value when the complex phase of the gauge coupling is 
large. As we emphasized in the introduction, the convergence proof of Refs, [y, [?J does not 
necessarily apply when the gauge cooling is employed; thus there is no contradiction even if 
the method leads to a wrong answer. Nevertheless, it is not yet clear what causes the failure 
for the gauge coupling with a large complex phase. To find the resolution of the problem we 
found in the present study, first we have to pin down what the real source of the failure is. For 
this, consideration on the basis of another approach to the functional integral with complex 
measure, the Lefschetz thimble [21-25|, might provide useful insight. See also Ref. [26| for 
suggestive observations. 
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